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CONVECTION EFFECTS ON RADIAL SEGREGATION AND 
CRYSTAL MELT INTERFACE IN VERTICAL BRIDGMAN 

GROWTH 


S. Tanveer* 

Mathematics Department 
Ohio State University 
Columbus, OH 43210 


ABSTRACT 

We analytically study the influence of convection caused by horizontal heat transfer 
through the sides of a vertical Bridgman apparatus. We consider the case when the heat 
transfer across the side walls is small so that the resulting interfacial deformation and fluid 
velocities are also small. This allows us to linearize the Navier-Stokes equations and express 
the interfacial conditions about a planar interface through a Taylor expansion. Using a no 
tangential stress conditions on the side walls, asymptotic expressions for both the interfacial 
slope and radial segregation at the crystal-melt interface are obtained in closed form in 
the limit of large thermal Rayleigh number. It is suggested that these can be reduced by 
appropriately controlling a specific heat transfer property at the edge of the insulation zone 
in the solid side. 
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1 Introduction 

In a vertical Bridgman apparatus, a cylindrical container (Fig. 1) containing a melt of a 
binary alloy mixture is translated downwards from a hot to a cold zone so as to cause so- 
lidification. The uniformity of composition of the resulting crystal and the relative absence 
of crystal defects are desirable features for technological applications. Ideally, these can be 
achieved if convection is eliminated and the crystal-melt interface is planar. However, in 
practice, this is difficult to ensure. Considerations of constitutional supercooling (morpho- 
logical instability) and the need to avoid transient effects due to container ends require a 
relatively large temperature gradient. On the other hand, a completely one dimensional 
imposed temperature gradient that would occur (Tiller et al 1 ) when the cylinder sides are 
insulated would require unrealistic large temperature differences between the two cylinder 
ends given that the length of the cylinder has to be large enough to avoid transient effects. 
Thus, one is forced to a configuration where significant heat flux occurs through the cylinder 
sides. However, when this happens, fluid next to the cylinder ends is hotter than the fluid 
at the center leading to convection for any value of the solutal or thermal Rayleigh num- 
ber. The review by Brown 2 discusses the Bridgman problem in great detail. Other review 
papers 3-5 deal with various aspects of directional solidification in general. 

There are many papers in the literature that address the problem of onset of convection 
in finite and infinite geometries (see references 6-15 and references there in). In these cases, 
the equations allow for a basic quiescent state (no fluid flow) that is stable upto a certain 
critical Rayleigh number. Such convection can be termed natural convection, as opposed 
to the induced convection caused by heat transfer through the side walls of a Bridgman 
apparatus. Since this paper concerns situations which in the absence of horizontal heat 
transfer is a thermally stable configuration, natural convection is not relevant in this context 
unlike that of an otherwise thermally or solutally unstable arrangement where possibility 
exists of a ’’resonant” response leading to vigorous convection. 

There are many papers dealing with the Bridgman problem itself (see Brown, 1988 for 
references). Nonetheless, the first set of fully consistent calculations that account for con- 
vection in the fluid, and its coupling with thermo-solutal field, heat diffusion in the solid 
and a non-planar interface appear to be due to Chang k Brown 16 and Adornato k Brown 17 . 
By assuming a quasi-steady growth when the transient scale is small, they calculate the 
nonlinear steady states numerically. Later on, these calculations were further extended by 
Kim k Brown 18 by including effects of heat transfer through an ampoule that surrounds the 
cylinder in an experiment. Unsteady transient effects have also been studied by Brown k 
Kim 19 and these calculations are in good agreement with experiment. These calculations 
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suggest among other things that the classic one-dimensional modelling through the Scheil 
equation 20 is an oversimplification since it assumes complete and through mixing all over 
the fluid. It is known for instance that between the diffusive limit, where convection can be 
totally neglected and the thoroughly mixed limit for which Scheil’s equation applies, there is 
an imperfect mixing zone at moderately large Rayleigh number. In this range, radial segre- 
gation and the interface deformation are larger than at very large Rayleigh number. Indeed, 
it was pointed out that reduction of gravity by a factor of 10 3 - 10 4 compared to earth 
may be detrimental to growing a crystal of uniform composition. Despite such advances in 
theoretical understanding, it is difficult from numerical results alone to get a global under- 
standing of the parameter dependences since there are so many of them (See Table 1). The 
trends in a certain subset of parameter space need not reflect the trend in other ranges of the 
parameter space. Thus, there is need for analytical results, which is likely to be of limited 
validity; nonetheless it can be complimentary to numerical calculations. 

The only analytical work that we are aware of that is relevant to convection in the 
Bridgman apparatus is due to Brattkus k Davis 21 . They analyzed a two dimensional model 
where the vertical dimension is far larger than the horizontal dimension and the heat flux 
through the side walls is assumed small. Brattkus k Davis 21 specifically concluded there 
was no necessary relation between radial segregation and mterfacial shape, a hypothesis put 
forward by Coriell k Sererka 22 based on diffusion alone. The numerical results of Brown and 
his coworkers, on the other hand, suggests a strong correlation between the two. However, 
the Brattkus- Davis analysis ignores the insulation zone shown in Figure 1. This insulation 
zone length is known to be an important parameter from prior numerical work (see Chang 
k Brown 16 , for instance). 

Here we consider steady state 2 inside a vertical cylinder (Figure 1) between 2 = 0 and 
2 = 2 2 , where 2 2 is assumed a constant. As in Brattkus k Davis 21 , the horizontal heat transfer 
is assumed small enough so that both the fluid velocities and the interfacial deformation are 
small. Prior numerical work (Chang k Brown 16 for instance) show that the latter assumption 
is valid even for relatively intense convection. The fluid velocity on the other hand has to 
be small so that nonlinearities in Navier-Stokes equation can be ignored. This assumption is 
clearly unrealistic in many experimental situations; however, as discussed in section 7, there 
are reasons to believe that the results on radial segregation and interfacial deformation will 
hold in part of the nonlinear regime as well. 

The above assumptions allow us to linearize about a basic one dimensional state and 
Taylor expand the interfacial boundary conditions about the original planar boundary. By 

sin reality, this is a quasi-steady state (see Brown (1988)); however for purposes of this paper, we make 
no such distinction. 
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modifying the no-slip boundary condition at the cylinder side walls, we find a modal rep- 
resentation for each of the stream function, radial temperature and concentration gradient 
in the melt that decouples each modes and reduces the problem to a set of finite ordinary 
differential equations for each mode. By making use of general expressions for the perturbed 
temperature and concentration fields in the solid, we find that the interfacial conditions can 
be expressed in terms of effective boundary conditions on the melt variables at the original 
planar interface. With the effective boundary conditions on a planar interface, we solve the 
linearized melt equations and find explicit expressions for interfacial deformation and radial 
segregation in the asymptotic limit of large thermal Rayleigh number R 7. Among other 
things, we find that there is a boundary layer near the interface that scales as where 

each of the solutal, thermal and fluid velocity field change rapidly. The radial derivative of 
the temperature field and the radial segregation in the crystal, each scale as provided 

|ZZ C | << Rj. . The coefficient of R T 1 ^ 6 in each scale decreases exponentially with the dis- 
tance of the interface from the end of the insulation zone in the solid side. In a certain range 
of parameters, consistent with many experimental conditions, our results suggest that radial 
dependence of the interface shape and concentration will be roughly given by the Bessel func- 
tion Jo(Air), where Ai is the first positive root of *To(A n ) = 0 (i.e. zero of «/i). Further, we 
find that in this case, the Coriell-Sererka (Coriell & Sererka 22 ) hypothesis of proportionality 
between interface slope and radial segregation is approximately valid, though the constant 
of proportionality is different from what these authors find with diffusion only. In this case, 
we also point out a specific condition on heat transfer near the solid end of the insulation 
zone, which when satisfied, will result in minimal interfacial slope and radial segregation for 
large R T . 

2 Mathematical Model 

In a Bridgman apparatus (Fig. 1), a binary melt is contained in a cylindrical container 
of radius a that is translated downwards with constant velocity -U a z , where 2 is a 
unit vector along the axis of the cylinder that opposes gravity, as shown in Figure 1. The 
concentration of solute (one of the two components of the binary mixture) at the top of the 
cylinder z = i 2 is c = c 2 . The top of the cylinder is maintained at temperature T = T 2 , 
while the bottom of the cylinder is at temperature T = Ti, which is significantly smaller 
than the melting temperature To of a planar interface. The density of the melt is assumed 
to be p po at temperature To and concentration c 2. This will be the reference density. We 
scale all lengths by a, all time scales by a/U a , all mass scales by p 0 a 3 and all temperatures 
by T 2 — T q , leaving us with non-dimensional quantities of Table 1. Concentration will remain 
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unsealed, as it is already nondimensional. Given the fair amount of algebraic manipulation, 
we prefer to use this set of parameters as it keeps the equations looking simpler. Off course, 
in the common Engineering literature, it is common to express these quantities in terms of 
Peclet number Pe, Reynolds number Re, Prandtl number Pr, etc. For the benefit of the 
readers, we have prepared Table 2 that expresses our non-standard parameters in terms of 
the more well-known dimensionless numbers. In section 7, we also discuss our results using 
the more standard notation to benefit the reader who is more interested in the concrete 

results than in the analysis. 

Using the nondimensional variables in Table 1, the top end of the cylinder, corresponding 
to z = z 2 , is maintained at nondimensional temperature T = 1 and the bottom end is 
maintained at temperature 7 \ with T, < 0 . The temperature gradient is assumed 

strong enough to avoid constitutional supercooling. The heat is allowed to flow through 
the sides of a cylinder for z 2 > * > z n in the melt zone and 0 < z < */ m the 

crystallized zone. The heat transfer across the side walls is such that the fluid velocities are 
slow and the interface deformation from a planar interface is small. The precise limitations 
placed on the size of the heat transfer by this assumption will be examined later in section 
6. For the present, it suffices to assume that heat transfer is sufficiently small. We introduce 
a cylindrical coordinate system (r,z) where r is the radial direction. No azimuthal 
angular variable is necessary as the flow variables are assumed to be axi-symmetnc. The 


mathematical equations for the steady solution in the melt are 

v- V T = k V 2 T, (!) 

v-V c = DV 2 c, (2) 

V . v v = - Vp + {aT + (c - c 2 ) } g z + v V 2 v, (3) 

V • v = 0 , ( 4 ) 


where T denotes the nondimensional temperature, c the concentration of solute (one of 
the components) relative to the total (measured in molar fraction) and v is the relative 
melt velocity (see Table 1) and v denotes the nondimensional kinematic viscosity (i.e 
inverse Reynolds number based on U a ). In equations (3) and (4), we have invoked the usual 
Boussinesq approximation, where the density variation due to change in temperature and 
concentration from reference values are included only in the forcing term on the right hand 
side of (3). Here, a ( > 0) and /? are the non-dimensional coefficients of volumetric 

expansion due to increased T and c respectively. When the solute density is smaller than 
the alloy, /? > 0. However, there is no restriction on the sign of in our current analysis. 
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On the crystal side, since there is no fluid motion and so 


v = 


Z. 


Thus the equations for temperature and concentration fields are given by 

d T. 


dz 

dc 


= k, V 2 T s , 


dz 


a = D, V 2 


( 5 ) 


( 6 ) 


( 7 ) 


Now comes the boundary conditions. Denoting the r and 2 components of v by ( u ,u>), 
we take the boundary condition on fluid velocity components at the sidewalls r = 1 as 


u = 0 , 

dw 

~dr = °" 


( 8 ) 
( 9 ) 

Equation (9) is a no-stress condition. A no slip condition w = -1 would be more realistic; 
however, we could find no simple basis representation of the solution in this case. We suspect 
that aside from the changing the nature of a boundary layer near r = 1, the change of this 
boundary condition will be have no global effect on crystal shape and radial segregation, at 
least in the limit of large Rayleigh number. 

The condition of no mass flux through the side walls imply 

dc 


dr 


= 0 . 


Also, the condition of heat flux imply 

dT 


dr 


= 0 for /( 1) < 2 < ZlI , 


dT 

+ eq ' T 


t q 2 for z n < z < z 2 . 


On the solid side, at r = 1 , 


dT a 

dr 


= o 

dr 

= 0 for /(r) > z > z u 


while 


dT, 

dr 


+ e <J\,T, = e q 2 , for z\ > z > 0. 


( 10 ) 

(11a) 

(116) 

( 12 ) 

(13a) 

(136) 


Each of eqi , eq 2 , cq\, and cq 2 , will be taken as constants that characterize the heat 
transfer through the side walls. In the common engineering literature, tq x and eq lf are 


5 



Biot numbers related to heat transfer through the side walls in the melt and solid region 
respectively. For instance, in equation (16) of the Adornato & Brown 17 paper, the notation 
Bi(z) with a specific choice of a piecewise constant profile, is related to the constants c</i, 
and e< 7 i, while B*(*)M*) is related to eq 2 , and tq 2s , provided we ignore the ampoule in 
their model. The quantity e will be assumed small, while each of q i and q 2 assumed 0(1). 
c will be our perturbation parameter. The precise multiplicative decomposition of the Biot 
number e< 7 i into q\ and e is unimportant since in the final results, only the product e<7i 

appears. The same is true for tq 2 , e<7i, and tq 2t . 

At the cylinder top, i.e. 2 = z 2 , we have the nondimensional temperature 


T = 1, U 4 ) 

c = c 2 , (15) 

u = 0 , ( 16 ) 

w = —1, (1^) 

with c 2 being assumed a constant. Similarly at the bottom, z = 0 , we get 

T a = T u ( 18 ) 


where Tj < 0 and assumed smaller than the nondimensional melting temperature at the 
interface. Note that one can use more general boundary conditions than (19), however this 
appears to be most relevant since in the limit of an infinitely long cylinder, conservation of 
mass dictates that the concentration should be the same at z = ±00 . 

Now comes the boundary conditions on the solid melt interface 2 = f{r) . We will 
assume local thermodynamic equilibrium since the relaxation time of the departure from 
equilibrium can safely be assumed to be much smaller than the typical time scales in this 
problem. Thus, at the interface, we satisfy the melting condition given by the solidus- 
liquidus line with the incorporation of the Gibbs-Thompson effect (i.e. lowering of the 
melting temperature in the presence of curvature effects) 


T = T, — — mc + 


do | 


/" 


1 

+ - 


r 


(1 + ff 2 r (1 + r 2 ) 


,2x1/2 


}• 


( 20 ) 


Further, from the solidus liquidus line, we must have 

c s = k c. 


( 21 ) 
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The continuity of heat flow implies that 


T v ■ h — 


dT 

dn 


= -T.z-n - 


K s 


on 

dn 


+ S t v ■ h, 


( 22 ) 


where h denotes a unit vector normal to the melt-solid interface pointing towards the liquid, 
anc ^ dn denotes component of gradient in that direction, St is the nondimensional latent heat 
(Stefan number). From the conservation of solute concentration across the interface, we must 
have 

- - n dc — r, dc $ 

cv-n-D — = -c, z-n-D, . (23) 

Further, from integrating (4) over a small control volume centered about a point on the in- 
terface, it follows that normal velocity is continuous when the difference of densities between 
the melt and the crystal is neglected. Thus, 


v • h = — z * h. (24) 

Difference of densities between solid and liquid side of the interface can be accommodated 
by including an additional — i ps ~ p ) z ■ fi term on the right hand side of (24); however, here 
we will be primarily interested in situations where this is not important. Further, at the 
interface, if we assume a no-slip condition for the fluid flow relative to the solid, then 

u + (u> + 1)/' = 0. (25) 

While the no-slip condition (25) at the melt-solid interface seems to be widely used in 
the literature, we are unaware of any convincing physical argument this is superior to a 
no tangential stress boundary condition. It may be pointed out that unlike the case of 
a fluid next to a solid where good experimental evidence exists for the no-slip boundary 
condition, the liquid molecules in this case do not preserve their identity as it goes through a 
phase transition. Nonetheless, despite this uncertainty, we find that the asymptotic scalings 
and parameter dependence presented later in this paper remain the same with a no-stress 
boundary condition, only the scaling constants differ. 

Note that the mathematical model with equations (l)-(4), (6), (7) and boundary and 
interfacial conditions (8)-(25) contain twenty irreducible nondimensional parameters D, D s , 
k, k s , ag , fig, v , z 2 , zj, zu, eq x , eq Xs , eq 2 , eq 2s , c 2 , 7\, d 0 , m, k and S t . The relation of the 
less standard parameters in this list with the conventional Peclet number, Prandtl number, 
Reynolds number, etc. is given in Table 2. 


3 Steady state for zero epsilon 

Despite the complexity of the general equations, there exists well known simple solutions to 
the above equations, as determined originally by Tiller et al 1 . If e = 0 , a quiescent state 
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is a simple solution to the above equations in which the melt velocity v = - z and the 
melt-solid interface is planar, i.e. / = z 0 , a constant. The temperature and concentration 
fields in this case are denoted by a superscript 0 as they are the leading term of an expansion 
for small c . They are given by 


T° = Ti + [e - - lj e _ i( „_ wl _ 

(26) 

c° = Ci + (e ] - l) e _i( Z2 - Z0 ) _ J- 

(27) 

On the solid side, zq > z > 0 , 


= T, + - l) ei , _ 

(28) 

<« = - o _ ;• 

(29) 

where c, is the concentration value on the melt side of the interface and T, is the interfacial 
temperature. To determine the three constants T, , c, and z 0 , we use the boundary 
conditions (20), (22) and (23) which in this case simplifies to 

Ti = — m Ci , 

(30) 

, 1 - T< - T ; ~ — + s ti 

e --(*2-zo) _ 1 e K ’ Z ° - 1 

(31) 

02 - Ci I. - C 2 - k Ci 

c ' e -i(v- z o) _ i ' e Tk z ° - l’ 

(32) 

In this case, the pressure p = p°(z) is hydrostatic and satisfies 


= a gT° + (3g(c°-c 2 ). 
dz 

(33) 

Note that all other equations and boundary conditions are trivially satisfied by the solutions 

(26)-(29) and (33) provided they satisfy (30)-(32). 

An important limiting case (valid in many experiments) is that k/z 2 

>> 1, K, jz 2 » 1. 

This simplifies expression for T° and T a ° as 


T° ~ —me, T (1 + me,), 

z 2 - z 0 

(34) 

T° ~ — mq — ( T\ + mc(), 

* Z 0 

(35) 
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each of which is a linear function of 2 . In addition, when (z 2 — zq) » D, S t << k/ z 2 , 

S t « k s /z 2 and D s « z 0 , (conditions that appear to be valid in some experiments), 
the matching conditions imply that the interfacial concentration 

1 

c ‘ ~ (36) 


Zq ~ 


Zl < 1 + 


*0 + T<») l" 1 
K s(~J LC 2 ~T \ ) J 


(37) 


Notice that the quantity is the ratio of the temperature difference between the 

top and the interface and the temperature difference between the interface and the bottom. 
Clearly, by controlling the location of the insulation zone, z 0 can be made to lie between 
z = Zj and z = 27 /, as will be assumed here. 


Until we get to a discussion of the concrete formulae ( 212 ) and (213) for interfacial shape 
and radial segregation in sections 6 and 7, the only simplification that will be used is that 
each of T° and T° is a linear linear functions of £ as in ( 34 ) and (35), since k/z 2 and k s /z 2 
are assumed large. Otherwise, the analysis will proceed with the assumption that each of 
c°, T°, c° and T s ° are known from (26)- (32) . 


4 Perturbed Steady state for nonzero e 

Now consider a small nonzero c . In this case, the dependent variables can no longer just 
depend on z . As is well known, the presence of a radial thermal gradient means that the 


quiescent state is no longer a steady 

state solution 

to the problem. 

We express solutions as 

a perturbation expansion in powers < 

e : 




T 

it 

o 

+ 

t T 1 


(38) 

c 

+ 

o 

V 

II 

c c 1 


(39) 


u = t u 

+ 


(40) 

XV 

= -1 + 

e w 1 

+ -1 

(41) 

p 

= P° + 

ep' 


(42) 

T, 

II 

+ 

tT] 

+ 

(43) 

C s 

= c? + 


+ -‘5 

(44) 

f = 

= Zo + e 

f + ... . 

(45) 
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Substituting these into (1-4), (6) and (7), and equating the e terms in the resulting equations, 


du l 

dz 


dv u 1 
dz 


On the solid side, we have 


c* 

> 

54 

II 

© 

+ 

CO ^ 

(46) 

— + til 1 c 0 ' = D V 2 c 1 , 
dz 

(47) 

= - a -f + r v» »' - V, 

or r 

(48) 

— + ( a T 1 + /? c 1 ) g + i/ V 2 w 1 , 
dz 

(49) 

- zr{ru l ) + ^-ta 1 = 0. 

r dr v dz 

(50) 

_ r yZ rp\ 

— K s v 1 3 , 
dz 

(51) 

JA = D.V’c]. 

(52) 


We note that in this coordinate system, the operator V 2 is given by 


d 2 


1 d 


d 2 


V 2 - — — + - — + — . 

v ~ dr* r dr dz * 

The boundary conditions on r = 1 for z 0 < z < z 2 are 

u 1 = 0, 

d w 1 


dr 
d< "} 
dr 


= 0, 
= 0, 


dT_ 

dr 


= 0 for zq < z < zji, 


dT 1 

dr 


+ <71 T° = 92 for zil < z < z 2 . 


On r = 1 for 0 < z < z 0 , from (12)-(13), we have the boundary conditions 


dT]_ 

dr 


dc 1 

-t = 0 , 

dr 

dT 1 

— - = 0 for z 0 > z > zj, 

dr 

= - T° + 92 , for Zj > z > 0. 


(53) 

(54) 

(55) 
(56a) 
(566) 

(57) 

(58a) 

(586) 
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The boundary conditions at the top z = z 2 are 



T 1 = 0, 

(59) 


c 1 = 0, 

(60) 


u 1 = 0, 

(61) 


w 1 = 0, 

(62) 

and at the bottom, z = 0 , 

II 

JO 

(63) 


<£ = o. 

(64) 

On the original interface z = z 0 , (20) implies 


T) + P if = T 1 + p T 0> = 

-m c 0 ' Z 1 - m c 1 + d 0 j/ 1 " + £ /*' } . 

(65) 

From (21)-(25), we get 



+ P c 

°,' = kc' + kp c°\ 

(66) 

—k T°" f x r ^ , T°" f l 

f K Tz K - T ‘ f K - dz ’ 

(67) 

-D c°" /* - D d A + (*_i) 

c> + c°' /■ ] = -n, cf /■ - D, 

(68) 


w l — 0, 

(69) 


u x = 0, 

(70) 

respectively. We can eliminate pressure from (48) and (49). Also, from (50), it follows that 


v (u ,w ) — curl A for some vector field A that can be chosen to be divergence free. For 
the axisymmetric flow under consideration A can be reduced to a scalar ‘stream function’ 
^ so that the velocity components of v x can be written as 


Then, from (47) and (48), we get 


u 


dtp 

dz' 


1 d(rtp ) 
r dr 


d_ 

dz 



xp 




Ip + g [a 


dT 1 

Ih 




( 71 ) 

(72) 


(73) 
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where the differential operator C is defined such that 

d 


C u 1 


dr 


f 1 d(r u) 1 
\r dr ) 


obtain 




d dc 1 
dz dr 


-O' 


(- 




(74) 

| and 

(47) with respect to 

r , we 

d 2 \ 

1 T. 

(75) 

dz 2 j 

' dr 

d 2 \ 

1 dc 1 

(76) 

dz 2 j 

1 dr 

, to r , 

we get 



d 

dT} 

(r ^ d " \ 

dT} 

dz 

3 — k 

dr 

+ dz 2 j 

dr ’ 

d 

dc\ „ 

, 92 ^ 


'dz 

— = 
dr 

\ C + dz 2 ; 

I dr' 


The boundary conditions (53) and (54) imply that without any loss of generality, 

(M) = 

£ VKM) = 

while the boundary conditions (61), (62), (69) and (70) imply that 


(77) 

(78) 

(79) 

(80) 




t/i(r,z 2 ) = 

=0, 

(81) 



9 *l> , \ 

& (r ’ Z2) 

= 0, 

(82) 



(»’, 2o) 

= 0, 

(83) 



9 i> , \ 

(r ’* o) 

= 0. 

(84) 

Further on taking the tangential 

derivative (i.e. 

derivative with respect to 

r ) of each of the 

equations (65)-(68), we 

obtain 





ar? 

dr 

+ /''if = 

dj 1 + r o' 5 

dr 

(85) 

dT 1 
di ’ 

+ P'T- 

= -m c 0 ' j 

1/ — rn -7— + d 0 C , 

dr 

(86) 


dr 

+ r' c°: = k 

^ + k /■' C°', 
dr 

(87) 
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( 88 ) 


—k T°" f 1 ' — — 

f dz dr 


= ~K S Tf f 1 ' - K s 


d_ dT[ 
dz dr ’ 




— + C°' f 1 ' 
dr J 


= -D s cf f‘ 


D - Yz d 4- ^ 


In addition to the above boundary conditions, elementary considerations of smoothness of 
each of the variables T x , T] , c 1 , c\ and v x in the neighborhood of r = 0 together 
with consideration of axisymmetrical flow, leads us to conclude that as r -+ 0 , 


0 ~ fi (2) r, 

(90) 

dr 1 

dr ~ M*)r, 

(91) 

dr, 1 

(92) 

dc 1 

— ~ <4(2) r, 

(93) 

dc] 

~ < 5 (^) r, 

(94) 

.V 

y ~ constant r 

(95) 


for some functions t j through , their precise form being unimportant. 


5 Series representation of solution 


The problem (equations (73), (75)-(78)) in conjunction with (55)-(60), (63), (64), (79)-(89) 
and assumptions (90)-(95) completely determine each of the unknown functions 4~T X — c 1 

d l d f or ’ dr ’ 

dr^s i q^ c s an< ^ V’ as functions of r and z and / 1 as a function of r. 

We expand each of the variables 


dT x “ 2 

~ h w ° n {z) Jl(An r) ’ 


dc 1 a ^ 2 


= p hn ^ Jl ^ An r) ’ 

. ga ~ 2 

^ ~ — 2^ 777TT 0,1 (*) r )> 

^ n=l J \\ A n) 

dT] ~ 2 

dr “ ^ j;(A n ) ^ r *’ 

8c] a “ 2 

dr “ ^ 4 ” M J, * A " r) ’ 


(96) 

(97) 

(98) 

(99) 
( 100 ) 
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where A n is the n th 
representation since 


f v = E 


2 

•A'(An) 


d n "A(An 


( 101 ) 


positive zero of the Bessel function J A . Notice that this is a suitable 


C J,(A„ r) = -A* J j (A n r) 


and the boundary conditions (55), (57), (79) and (80) are automatically satisfied by each 
mode. The representation of ^ and in (96) and (99) may appear to contradict (56b) 
and (58b). This is not the case since the convergence of the series in (96) and (99) for 
z > zji and z < z; is only in the mean and not pointwise at r = 1. In this range of z, 
there is a slow decay of a n and a, n as n -+ oo, which implies that one cannot calculate 
the derivative of the series (96) and (99) term by term. Thus, it is not possible to plug the 
series expressions directly into the differential equations and obtain the correct equations 
for a n and a sn . Instead, we multiply each of (73), (75)-(78), (85)-(89) by r J, (A n r) and 
integrate with respect to r from 0 to 1. On integrating by parts, and using (55)-(58), (79), 
(80), (90)-(95) and the following relations: 

a " <2) = ik l dr r r) V (r,2) ’ 

« = dk l dr r J,{Xn r) * (r,z) ’ 

<=«(*) = g a JiM /» 

•-W - lk) l 0 irrJ '^ T ^^ z) ' 

< 2 > = dk) /o dT r MK r) (r,2) ’ 

d n = 1 T t dr r Ji( A n r) /(r), 

we obtain three ordinary differential equations for each n on the melt side and two ordinary 
differential equations on the solid side. These equations can be written as 


C 2 b n = X 2 n R c h(z) c n for z 0 < z < z 2 , (102) 

£ 3 Cn = -a n - b n for z 0 < z < z 2 , (103) 

£ lm a n = -X 2 n Rt c n T A n [ -<?i T° + q 2 ] for z u < z < z 2 , (104a) 

£ lm a n = -A* R t c„ for z 0 < z < z n , (1046) 
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A, b sn = 0 , for 0 < z < z 0 , 


(105) 

A, dsn = A n [ —q it T° + q 2s ] for 0 < 

z < zi 

(106a) 

A, a sn = 0 for Zi < z < z 0 


(1066) 

where 



N 

1 

N 

T 

II 

"bT 


(107) 

II 

* 4 
p 


(108) 

R _ 9p ( d ~ C 2 \ 

Dv \1 — e _ l 22-2 °)/ D / 


(109) 

and the differential operators 



£l " = d? + K Tz ~ A ” ’ A ‘ = + 

1 ^ >2 

(110) 

Cl = D J? + Tz - x * D ' 


(111) 


d 

dz' 

(112) 

+ E - A » D - 


(113) 


Note that since k / z 2 » 1 , (34) holds and therefore T 0 ' is a constant. So the thermal 
Rayleigh number Ft? defined in (108) is also a constant (Note the definition in terms of 
dimensional variables as well in Table 2). 

From (59)-(60), (81) and (82), we find that at the cylinder top, 


a n (z 2 ) = 0, (114) 

K (*a) =0, (115) 

C n (Z 2 ) = 0, (116) 

«n {**) = 0. (117) 


The boundary conditions (63), (64) at the cylinder bottom imply 


°»n(0) = 0, (118) 

Ki (0) = 0. (119) 

At 2 = z Q , we find from (85)-(89) that 

a sn(zo) + d n T° {z 0 ) = a n (z 0 ) + d n T^izo), (120) 
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a n {z 0 ) + 4 T°'{z 0 ) = ~m c°'(z 0 ) d n - m - b n {z 0 ) - 4 X 2 n 4, (121) 

b sn (z 0 ) + d n ^ c°'(z 0 ) = k b n (z 0 ) + kd n £ c°'(z 0 ), ( 122 ) 

— /c < (z 0 ) - K d n T°"{z 0 ) = -K, a sn (z 0 ) - 4 4°"(z 0 ), (123) 

-Db' (z 0 ) - Dd n -c°'\zo)+(k-l)[b n {zo)+c 0 '(z 0 )^d n ] = -D s b' sn ( z 0 )-D s d n -c° a (z 0 ). (124) 
From (83) and (84), it follows that 

Cn(z 0 ) = o, (125) 

4 (zo) = 0. (126) 


Before proceeding further, for purposes of reducing algebra, we simplify the equations 
further by assuming that 

(a) K, k/z 2 , Ks and k s Iz 2 » 1 
(b) 4=0 

With the assumption (a) as above, we replace equations (104a, b) and (106a, b) by simpler 
equations 

£i a n = -A* R T c n + A n [ -qi T° + q 2 ] for z// < z < z 2 , (127a) 

Ci a n = -A* R t Cn for z 0 < z < z /7 , (1276) 

Ci a sn = An [ -qi. 4° + q 2 ,] for 0 < z < z; , (128a) 

Ci a an = 0 for zj < z < z 0 , (1286) 


where 


Cl 



(129) 


The interfacial matching conditions (120)-(124) couple the concentration and temperature 
fields in the solid side to the melt variables. Using (105), (118), (119) and (128), expression 
for a s1l and b sn can be obtained in terms of two arbitrary constants. By eliminating these 
arbitrary constants between a sn (z 0 ), a' sn (z 0 ) and 6 SJl (z 0 ), 6' an (z 0 ), the matching conditions 
( 1 20)-( 1 24) can be written as two effective boundary conditions at z = z 0 on the melt 
variables. As shown in Appendix A, these boundary conditions are 


— — a»i(zo) T 4 fln(-o) 4“ fcbn(zo) 

A n 


A* 


— ^n(^0 z l) 
c 7 


b' n (zo) 4" m \b n {z o) 4" ^ 2 ° n (^o) 0? 


( 130 ) 

( 131 ) 


where the effective interfacial parameters fa, fa, rn,, m, and fa are defined as 

_ k,(T? + mc»'(z„))(l + e-“"'") 

11 ~ a(mc»'(ao) + r°')(l )’ 


( 132 ) 
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From (102), (103) and (127), we get an eight order linear differential equation for c„ : 

[£i£ 2 £ 3 + \ n R c C\h(z) — \^RtC 2 c n = £ 2 [— A n [— q^T° -f </ 2 ]j for zn < z < z 2 , 

(143a) 

[£i£ 2 £ 3 + A n R c C\h(z) — A^/? 7 ’£ 2 J c n = 0 for zo < z < zn- (1436) 

Alternatively, instead of one equation for c n , we can write the following equations for c n and 
a n 

[^£2 — D\ n Rx + \^R c h(z) c n + a' n = — £>A n [— q^T 0 + </ 2 ] for zjj < 2 < z 2 , (144a) 

[£ 3 £ 2 - D\ 2 n R T + \ 2 n R c h(z) c n + a' n = 0 for z 0 < z < z f/ , 

£1 a « = —^ n R-TC n + A n [—qiT° + q 2 ] for zjj < z < z 2 , 
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(1446) 

(145a) 


(1456) 


C\ a n = —X^RtCu for z 0 < z < z ii ■ 

The boundary conditions are (125), (126), (130) and (131) at z = z 0 and ( 1 14)-(1 17) at 
2 = z . 2 . In addition, we require that a n and its first derivative, b n and its first derivative 
and c„ and all its first three derivatives are all continuous at z = z u . From (103), this is 
equivalent to requiring that a n and its first derivative, c,, and its first five derivatives are 
all continuous at z = z//. These are enforced on a suitable representation of the general 
solution, as will be discussed shortly. 

Since (143a) is an eighth order non homogeneous linear equation, we can express the 
general solution for zjj < z < z<i in the form 

C. = + t^ C > F ‘(z) + ErtW, ( 14 «) 

j = i i= 3 

where c np (z) is a particular solution to (143a), F/ are eight independent solutions to the 
homogeneous equation (143b) and 

7 = a^ 3 4 /6 - ( 147 ) 


Factors containing powers of 7 in (146) can be absorbed as part of arbitrary constant C } 
by suitable redefinition. We choose not to do so for we want to explicitly show the R T scaling 
of the eventual answer. With our choice of Cj, it will turn out Cj = 0(1) as Rj 00 • 
Questions of determination of functions F/ and c„ p will be set aside for the moment. Once 
c„ is obtained, expression for a n and 6 n be found, at least in principle from (144), (145) and 
(103). We write these symbolically in the form 

an=a np (z) + 't-f'C i F;(z), (148) 

j = l 


2 \2 O 8 

K = K, + E C’CyJjW + -br E OiF }(z ) , 

j=l > 3 = 3 


(149) 


where F] are related to while particular solutions a„ p and b„ p are determined in 
terms of c„ p . 

For z 0 < z < zjj , i.e. in the insulated zone, we write 


c . = 'Er 7 c i Fj + Y.T'Ci F i 

j = i i - 3 


(150) 


a„ = Y,r'CjFj 

i= 1 


(151) 
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(152) 


2 \2d 8 

6„ = £ 7 -vv?+^ £<?,*? 

J=1 ' j=3 

In appendix B, we derive asymptotic representation for each of F c , F a F b c a 
and 6 np for large that is uniformly valid for all A n , provided \R C \ « R T . We also 
derive simpler expressions for the special case R T » 1 with 7 s >> A 2 n R c and this 

latter expressions will be quoted here as these are the only ones needed in obtaining concrete 
results for |/? c | << R^ 6 . 

In the special case, A 2 |/? c | << 7 s , the leading order asymptotic behavior is given by: 



c A„[- 9l r" + , 2 ) 

[A^r + AfJ y{\ 2 n R T + xiY' 

(153) 


F, C ~ — e P i ( z ~ z °) fc ^ El e P 2 (z-zi) 
D 2 D 

(154) 


F- ~ c 7i(*-*o) f or 3 , 4,5 , 

(155) 


FJ ~ e^- 2 *) for j = 6, 7, 8 , 

(156) 

where pi and P 2 are given by 



2D*\l4D> +X - 

(157) 

Also, 7 j for j =3 to 8 

are the six independent roots of 



L Aii\ L Aij] - ^Rt = 0 , 

(158) 

where 

<N £ 
1 

CN 

S* 

II 

s 


The roots of (158) are 

I'M = (y 2 - A If + y -(y 2 - XI). 

labelled such that for Rt » A*, 



7j ~ lUj 

(159) 

where 

UJ 3 — — e 1 = —cog 

(160) 


OJ 4 — — 1 — —uJy 

(161) 


cv 5 = -e ~ iv/3 = -w 6 

(162) 
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a np is general is a complicated expression, given in the appendix. We will not need the 
general expression in what follows except to not that in the special case when A* « R T , 
the general expression simplifies such that 


a 


up 


A|Hir_+d A t 5 1<7 iT 0 '/?t A lRch{z)_ 

XlRr + Xt v[\\Rt + F n] 2 L 2 [—D *] 


,+ 


np 


npT f.j-n - ’I 


(163) 


where 

L 2 [y] = D{y 2 - A^) + y- 


For 7 5 >> Xl\Rc\ (A„ unrestricted otherwise), 



p° ^ e Pl(z - zo) F£ ~ e P2 * z 

-z 2 ) 

(164) 

*7 

^ _ "7 for j = 

l] ~ K 

3,4,5 , 

(165) 

F j 

~ "7 e Tf>(*-*a) f or j 

-v? _ a 2 

1 j n 

6,7,8 . 

(166) 

The asymptotic form of b np 

is simply obtained by using 




bj\p z — ^np ^3[^np] 


(167) 


For 7 5 >> A^|f? c |, The independent homogeneous solutions for b n have the asymptotic 

behavior 

pb ~ _ e Pi(z-z 0 ) , pb ~ _ e P2 (z-z 2 ) , (168) 

pt> ~ _piA z ~ z o) l DL ^} - h(z) + .. for j = 3,4,5 , (169) 

J 7^2 rel 

p b ~ _ e -rj {z-z 2 ) 2Rk^lA h(z) + .. for j =6,7,8 , (170) 

3 

where 1 

These asymptotic results for a special case will be used later in obtaining concrete ex- 
pressions for quantities of interest. It is to be noted that the asymptotic expressions (155), 

(156), (165), (166), (169), (170) become invalid when \X 2 n R c \ is the same order as 7 5 or larger. 

This can occur even for R c « R-T provided A n is large enough. In Appendix B, section 
2, we derive more general expressions for the solution of the homogeneous equation will be 
uniformly valid for all A n provided R T « Rc The expressions (155), (156), (165), (166), 
(169) and (170) will have to be replaced by these expressions. The remaining asymptotic 
relations (153), (154), (163), (164), (167) and (168), however, remain valid in this range of 

parameters. 
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Despite differing asymptotics in different regimes, at this stage we prefer to think of 
relations (146)-(152) as exact. The only aspects of the asymptotic results of Appendix B in 
differing regimes that will be used in this part of our analysis is that in all cases, 

(1) Each of Fj, Fj and Ff have a common exponential part multiplied by an algebraic 

dependent term that depend on superscript a, b or c. For instance, when | | << 7 5 , 

Wl = p '( z ~ z °)’ ^ = P 2 (z~z 2 ), Wj = 7j (z-zo) for j =3,4,5 W, = 7 j (z-z 2 ) for j = 6,7,8 

(171) 

(ii) The variation in z of e w > for j = 3, 4. .8 is much larger than the variation of the algebraic 
prefactors when R T » 1, while for A n >> 1, the variation of the exponential part e W} 
for j = 1,2 dominate the variation of the prefactors. 

(ii*) For j = 1,3,4, 5, Re Wj > 0, while for j = 2, 6, 7, 8, Re Wj < 0. 

0 V ) The choice of the constants Cj and Cj and the prefactors is such that Wj(zo) = 0 for 
j = 1,3, 4, 5, while Wj(z 2 ) = 0 for j = 2,6, 7,8. 

As a consequence of the above properties, it is clear that for large R T , we have each of 
F j( z o)> Fj( z o) and Fj{z Q ) will be exponentially small in R T for j = 6,7,8, while FHz 2 ), 
Fj (z 2 ) and F '■ (z 2 ) are exponentially small in Rj for j = 3,4,5. Again for A n >> 1, each of 
F 2 ( z o), F>(Z 0 ) , F£(zo), Ff(z 2 ), Fj(z 2 ) and Fj(z 2 ) is exponentially small in A n . These facts 
will be needed later in simplifying matrix equations. 

1 he continuity of a n and its first derivative, c n and its first five derivatives can be written 
as the following matrix relation 


GX = R, (172) 

where G is a 8 x 8 matrix, X and R are column vectors with 8 entries. These are defined 
such that for j = 1,2, 



Xi=T\6 i -C i )e w *’"\ 

(173) 


S’-'t 



3 ~ dz«-') z=ZlJ an "’ 

(174) 


F i f» r k - 1,2 , 

(175) 

and 

for k = 3, 4, ..8. 

(176) 

For j = 3, ...8, 




Xi = {Ci - 

R — 7 8 ~j ^ (j 1 c 

3 7 dzU-i) 


G kJ = 


d< k ~') 

dzG~ l ) 


F“ for k = 1 , 2 , 


(177) 

(178) 

(179) 
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and 


(180) 


G kj =7- ( ‘- 3 »e- W '(->^)_ i /J for k = 3,4,..8. 

Note that in arriving at (172)-(180), we had to scale the original continuity equations with 
appropriate powers of 7. Further, we scaled the unknown X with e w >( z "l These steps are 
necessary to ensure that the limiting G as Rt — * 00 0 e - 7 00 ) 1S 11011 s ’ n S u ^ ar 

free of transcendental terms in Rt- 
Denoting 

H = G~\ ( 181 ) 

it is clear that we can then write for j = 1,2, 

6 j = C, + E (182) 


k= 1 


while for j = 3, ..8, we have 


Cj = Cj + E 


(183) 


Jt=l 


We satisfy (125), (126), (130) and (131) at 2 = z 0 and (114)-(117) at 0 - z 2 , which on using 
(182) and (183), is equivalent to the matrix relation 


MZ = 5, 


(184) 


where M is a 8 x 8 matrix, while Z and S are each column vectors with eight entries. These 
are defined such that 

Z = (Cj, C2, C3, C4, Cs, Cq, C7, Cs) ) ( 185 ) 

s = ( 0 , 0, P3e~ Xn{zo ~ Zl \ 0, 0, -7 a np (z 2 ), -7 5c « P ( z 2), -7 4 <4 p (z2), - 7 K p ( z i)) 


+ (Pl,/ , 2,P3,T4,Ps,/ , 6,F7,P 8 ) T 

where for 1 < / < 4, 

8 


p = - 7 m, E - E E R ^ R i 

j= i fc=6 i=i 

and for 5 < / < 8, 

P, = > E + E E H t ,,R r 

j= i fc=3 i =1 

The elements of the matrix M can be written as 

Mij = 7" 2 P/(zo), forj = l-,2 and M,j = Ff(zo) for j = 3, 4.. ,8, 


(186) 


(187) 


(188) 


(189) 
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M 2 j = 7~ 3 Ff(zo), for j = 1,2 and M 2j = ^Ff (z 0 ) forj = 3,4.,8 (190) 

M 3 j = 7' 1 \K'F- + A/? + fl 2 F) J 2o for j = 1,2, (191) 

^ 3 ,, = [ 7 ~ 1 a; 1 f;' + 7 -ViF; + ^a^ c 7 - 6 f ;] 20 for j = 3,4.., 8 , (192) 

= [a; 1 Ff + mi F- + m 2 F/J ^ for j = 1 , 2, (1 93) 

M 4,i = [Ff + m 1 Fj’] zo + m 2 F“(z 0 ) for j = 3, 4.., 8 , (194) 

M ^ = F-{z 2 ) for j = 1,2, .. 8 , (195) 

M e,j = 7 _2 F/(z 2 ), for j = 1,2 and M 6J = F^za) for j = 3, 4.., 8, (196) 

M 7j = 7 ~ 3 F/( 2 2 ),- for j = 1,2 and M 7J = 7 _1 Ff (z 2 ) for j = 3, 4.., 8, (197) 

~ Fj(z 2 ), forj = l,2 and M 8 j = A^R c7 _5 F] , (z 2 ) for j = 3, 4 .., 8. (198) 


Note that in arriving at the specific form of matrix M , which is nonsingular in the asymptotic 
limit of Rf > 00 , we needed to multiply each of the original boundary conditions by 
appropriate scale factors involving 7 . Once solution Z = M -1 S is found to (184), the 
quantities a n (z 0 ) and b n (z 0 ) can be written in terms of Z as: 


a n(z 0 ) = 7 -1 ^{z 0 )Z x + F 2 (20) | 


+ 7" 1 E F*{z 0 )Z } + 7- 1 E F°(z 0) Zi + e~ w *") E H hk R k ) , 


bn(z 0 ) = 7 1 j/?(2 0 )Z| + F%(z 0 ) [z 2 + ie - W2(zn) E H 2ik R k ^j | 

A 2 /? 5 A 2 /? 8 / 8 \ 

+-^r E fiMZj + -# E d(*o 2> + E «>,***) • (200) 


7 2=3 7 ” “ 6 ' V £S 7 

In terms of 0 ^( 20 ) and 6 n (^o)> as shown in Appendix A, the interfacial deformation coefficient 

a n (zo) + ^jfKizo) 

rf H -t \ 

'-*'71 


T« + mcO'(z 0 ) ’ 


and the radial segregation coefficient 


M*o) = (&c°>o) - C°'(2 0 )) + kb n (z 0 ). 


(201) 


( 202 ) 


Recall that the function e w > reflects the common exponential part of the growth for 
the functions F“ , Fj and FJ . As discussed earlier, and shown explicitly in appendix B, for 
R-t >> |F C |, with R t >> 1, e w J for j = 3,4,5 decreases exponentially at a large rate with 
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increasing z from a value of 1 at z = z 0 . This is uniformly true for all A n . The decrease is at 
least like Further, e w > for j =6,7,8 decreases exponentially with decreasing z 

from a value of 1 at z = z 2 at a rate that is at least 0{e'^ (z2 ~ z) ). Thus, the elements M ktj 
for 1 < k < 4, 6 < j < 8 and for 5 < k < 8, 3 < j < 5 are which is 

exponentially small in R T . Further, if A n » 1, it is clear from (157) that Pl A n and 

V2 ~ A n and so M ka for 1 < k < 4 and A# w for 5 < fc < 8 are each 0 (e~^- z ^ which is 
transcendentally small in A n . Thus, if A n or z 2 - z 0 is of the order of some positive power of 
R t , these terms will also be transcendentally small in R T . Further, using the growing and 
decreasing properties of e w G it is easy to see that the the elements P, defined in (187), (188) 

satisfy 

p l = 0 f e -M*//- 2 o) i e _ 5'’ , ( 2;/ -*°),e - *( Z2_; ' // ),e - ' Xn ( 22 " J!// )) • (203) 

Further, the term \~ x appearing in the third element of S in (186) is transcen- 

dentally small in A n for large A n . Noticing the special structure of the reduced matrix M, 
obtained by ignoring transcendentally small elements in the limit of large A n and R T , >t is 
clear that the solution Z to (184) has components 


Zj = O fe- A ' ,(2//_2o ^,e _ 2 ' v(2n-zo) ,e _An(2o_2/) ) for j = 1,3, 4, 5, 
Further, Z 2 = 0(1) and each for j = 6,7,8 are at best 0(1) and so 


(204) 


£>>'*(»> R,J = O , (205) 

e w '<“) = O (e-J (= "“ TO) ) for j = 6,7,8. (206) 

Thus, from (199)-(206), it follows that the effect of the insulation zone is to exponentially 
quench the small scale components of the radial segregation and interfacial deformation, 
which correspond to b sn (z 0 ) and d n for large A n . Most of the contribution in the summations 

(96)- (101) therefore comes from terms where A n = 0(1)- 

In that case, the matrices G and M simplify further in the asymptotic limit R T » T 
As shown in appendix C, (184) can be solved in closed form in this asymptotic limit resulting 

in (Z 6 , Z 7 , Zg) = o(l), 


(Z 3 ,Z 4 ,Z 5 ) 


1 .y/3 1 .V3 N 

2 + , T ’ ’“2 6 , 




-A„(2:o-z/) 


(207) 


-A n (2 0 --;) 


(Z,,Z 2 ) 


2m 2 /j 3 e 

(p, + mj - m2) - e6 ,1_PJ M ;!2_2o) (P2 + m > - m a)) 


(>• 


p P \{ z 2 


ro) ) , (208) 
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K(zo) ~ 27 1 0 3 e~ Xn( ‘ Zo ~ z ^. 


where the effective parameters /? 3 , mi and m 2 are as defined in ( 134)-( 136). From (199) and 
(200) and the above asymptotic results, 

«»(*,) ~ - 2 (, + m 2 [l-e<^><*-»>] | 

L [Pi + m, - m 2 - e (p>-P2K^- z o)(p 2 + mj _ m2 )]J 

M*0 ~ 2 7 -'/? 3 e- A "«'"-"'^ m 2 (l - ,.<»-»)(«■-■.)] 

[pi + m, - m 2 - e(P‘-P 2 H‘' 2 -‘'o)(p 2 + mi - m 2 )] 1 > 

From (201), we find 

, o/a -1 e~ Xn ( z °~ Zl '> [" m 2 (l — m^) (l — e( pi_P 2 H 2 2- 5 o)\ 

d » ~ 2^37 tt™?- 577^-17 1 + 7 : 7^ I F2111 

(i + me (20)) (pi -f mi — m 2 — (p 2 -f rri\ — m 2 )e( p i _p 2)( 2 2-zo)) ' > 

From (202), (210) and (211), we can obtain asymptotic expression for 6 4n (2 0 ) as well. From 
(4o), (101) and recalling that J((A„) = Jo(A n ), we get interfacial slope 


d n ~ 


2^37 1 


g — ^n{zQ Zj ) 

( T °’ + mc°'(z 0 )) 


(209) 


( 210 ) 


(211) 


MKr) 4<yj 3 _ An 

n=l \l /3 Jo(K) Rt 6 (T 0 ' + mc°'(z 0 )) 

m 2 (l — m|) |l — e t pi — p 2 ) J. 


— A n (^o — z 1 ) 


[ ' (Pi + mi - m 2 - (p 2 + mi - m 2 )e( p >- p 2)(*a-*o))J ‘ ( 212 ) 

Notice that in putting an upper limit 00 in the summation in (212), we have to go through 
an intermediate analysis where we replace the upper limit by N\, where Ni >> 1 but 

smaller than any power of R T since such a limitation is needed for the validity of (159) 
used in deriving (207) and (208). The contribution to the summation for A n even larger 
is transcendentally small, as argued earlier. Now since N x is far larger than unity and the 
series in (212) is convergent, the leading order asymptotics is indeed the same as with 
replaced by 00. From (44), (100) and (202), radial segregation in the solid at the interface is 


w s ( \ 
^r( r ^o) 


V) jWnr) 4e/3 3 (fcc°'(z 0 ) - fT(fo)) 
n=l An /3 J 0 (A n ) R]/ 6 (T 0 ' + mC°'(z 0 )) 


(Pi P 2 ) (^2 ^0 ) 


m 2 (l - m|) |l - e (pi-P 2 M 2 2-*o)j. 

(Pi + mi — m 2 — (p 2 + mi — m 2 )e( pi-p 2)( 2 2-zo)) 

+ E ^i(^n r ) 4 e0 3 m2kae~ Xn ^ z °~ Zl '> x |l — e ( p i-P2)(z2-z 0 ) j 

tlXl /3 MX n )0RT /6 ( Pl + mi - m 2 - (p 2 + m t - m 2 )e( p i- p 2)( 2 2-zo)) 


(213) 


Now, consider fluid velocities in different regions of the melt. By using (40), (41), (71), 
(72), (98), (108), (150), ( 154)-( 156), (182), (183), (185) and the asymptotic solution (207)- 
(208) and simple identities of Bessel functions, we can write the asymptotic expressions for 
radial and vertical fluid velocities for \z - z 0 1 = 0(R^ 1/G ) as 

A KRy 3 e0 3 e~ Xn ^ Zo ~ 2 i'> * .i/ 3 b i /6 < , 

U E— T^TTT— ~ E n T “ Az 2o) ^i(A„r) (214) 

n= 1 -An Jo[A n )l j — 3 
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(215) 


w + 1 


E 

71=1 


4 /6 efte' A ” (20 ^ 


Y Ar.e A " /3 ^ /6 ^ (2_ * o) Jo(A„r) 

X 2 n /3 J 0 (Xn)T°' U 


where 


.y/Z 

iV 3 = -!+* — 


.v/3 

N 4 = 2, JV 5 = -1 


This gives a boundary layer at * = z„ where the fluid velocity changes rapidly. Of course the 

numbers N, are such that at the interface a = So, each of these expressions reduce to « - 0 , 

u, = -l. Note that the expressions (214) and (215) imply that the velocities change rapidly 

within a boundary layer that scales as Aside from the factor eft, which generally 

depends on n and the heat transfer parameters, it is easily seen that within the boundary 

layer, the dimensionless horizontal velocity u scales approximately as T0 , e 

1 . * 1 Kfirp ~\i(z-zo) The 

the vertical nondimensional fluid velocity w + 1 scales approximately as -5#- 

form (214) and (215) suggests that there will be skinny convection cells within this boundary 

layer. More detailed discussions on the parameter dependence of velocities will be taken up 

in section 7. . , 

Now, let’s consider the neighborhood of 2 = z n , the end of the insulation zone in 

melt. We first consider z < z„. Since prior analysis shows that for j = 6,7,8, each of Zj 

and therefore Cj is at best 0(1), it follows that C 3 e w ^ is transcendentally small. Thus, 

from (153), (174), (178), ( 183 ) and (C15), it follows that for j = 6,7,8, 


e , \v 5 u n L * 1 v 

e Wj ( Z// )Cj ~ ~ H ji3 l c np {zn)) ~ 7 j,3 ^2 R T + \* 


A n [-<?i T°(zu) + <72] ^ (216) 


k - 1 


where to the leading order, it is known that % 3 = 0(1). Further, since 2, for , = 3,4 5 
turned out 0(1) at best, it follows that for j = 3,4,5, C,e^"> will be transcendentally 
small. Thus, in the expression (150), for the stream function coefficient, the leading order 
contribution comes only from j = 6,7,8 when z < *// with \z - z„\ small enough to be 
inside the boundary layer (actually an internal layer). It is clear from (150) that 


ye^ z - Zl,) H h3 ^c np (zn)) 


8 A n [- <jfi T°(z/j) + <72] pl ,(z-z u ) 
.3 ‘ 


j=6 


A IRt + A« 


Therefore, from (98) and (108), 


, ,, - 2 \ n KRrc\-<l\ T 0 (zi]) + 92) y. u 
~ ft(A„)(A ! „flr + K)T°' ” 


(217) 


n= 1 


Here, unlike what happens near z = z 0 , most of the contribution to the summation occurs 
when A n = 0(/?j 4 ) >> 1. In that case 7 j determined from (158) scales as R T and hence 
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from (71) and (72), we conclude that there is a layer that extends downwards from z = zjj 
of thickness R T ^ , where the radial velocity u scales as 


e*4 /4 [-g 1 r°(^ // ) + g2 ] 

rp 0 f 

and the vertical fluid velocity, w -f- 1 away from r — 1 scales as 

eK [~<h T°(zn) + <72] 


(218) 


(219) 


Near the corner, where each of 1 - r and 2 : - z n is 0(fl^ 1/4 ), the vertical velocity also scales 
as in (218), i.e. it is large as gravity effects become large. Similar arguments for z > zjj 
can be advanced to show that there is also a /?^ 1/4 boundary layer that extends upwards 
°f z ~ z " where horizontal velocities scale as in (218) and vertical fluid velocities scale as 
(219). Within this R T ' layer around z = zu, there are convection rolls that can be deduced 
from the r and 2 dependence in (217). We can put forward arguments near z = 2 2 as well 
to show that there is a similar boundary layer where convection occurs. 

Away from these boundary layers at 2 = z 0 , z = z n and 2 = z 2 , the horizontal and 
vertical fluid velocity components for 2 > z u scale as as can be deduced from (98), 
(108), (146) and (153). I his stays 0(1) as with increasing gravity effects. However, there is 
in this casual so a boundary layer near r = 1 of thickness 7^ 1/4 where the vertical velocity 
scales as and therefore intensifies with increasing gravity. 

For 2 < 2 //, but away from boundary layers at 2 = z u and 2 = z 0 , it follows from (150) 
and the behavior (155), (156) of Ff for j = 3, ...8 that these only gives transcendentally small 
contribution to c n . The only sizable contribution comes from possibly Ci and C 2 , which are 
determined in (208). Thus, it follows that in this region, the fluid velocity scales at best as 
T o'ftV« and is reduced by increasing gravity. Unlike the core region for 2 > z u there is no 
boundary layer in this case at the side walls near r = 1. 

In order for our leading order analysis to be self consistent, it is not necessary that scaled 
fluid velocities u and w + 1 be much smaller than unity. Indeed, we obtain the same leading 
order result if the advection term v- Vu were totally dropped in (3). In order for us to be 
able to linearize everywhere in the fluid field, it is necessary that the Reynolds number based 
on the largest fluid velocity Re } « 1. Based on the estimates of the velocities in the 
boundary layers at z = z 0 and 2 = 2 //, this would require that each of 


KRj 3 ( f 3 3 e ~ x '^ z °~ Zl ' 1 

1/T* 


« 1 


( 220 ) 


pT 0 ' 


« 1 


( 221 ) 
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For the complete validity of our analysis, e needs to be small enough so that both conditions 
(220) and (221) are met. In addition, we need that the interfacial slope, as given in (212), 
is small. The latter does not limit e particularly since for larger gravity levels, as on earth, 
this quantity is small for most experimental conditions without e being small. 

Equations (212) and (213) are general results on interfacial deformation and radial seg- 
regation that are valid whenever \R C \ « provided the constraints listed in table 

(2) are satisfied. However, this result is complicated by the dependence through effective 
parameters fc, m„ m 2 , which in turn are complicated functions of other parameters. It is 
therefore instructive to look at a special limit that is applicable to many experimental set 


up. 

6.1 Special case of D s jzo « 1 and A,(z 2 - Zo), A,*/. -Mo » 1 

In this case, as mentioned in section 2, « ~ £c 2 , and its derivatives negligible at 

* = c°'(z„) ~ -(l-l)g,c°>) ~ (1 - l)gr- Therefore, from (134)-(136), 


P 3 


mi 


K s 

K 


(i 


(9l,T s °(2/) ^ 

m(l — k)c 2 


m 2 ~ — ■ 


D D 2 (TO'-^(i-l))’ 
(3(\ - k)c 2 


"aD 2 (r°'-^(£- 1)7" 

Thus, in this case the interfacial slope expression (212) reduces to 

2D \/lD? £ 


dr 


^ 4J 1 (A n r)6 1 

“1 An /3 Jn(A T1 ) 


1 - 


-^n(zo-Zl) 




( 222 ) 


where 


<5, = 




8 > = 


9i. 


j 1 0 / 


['h.T‘Ji-1) - fl2,r 

<? (i - fe)c 2 (i - =y) 

63 ~ «D 2 (T°'-(|- l)mc 2 /D)‘ 

The expression (213) for radial segregation in the crystal at the interface becomes 


(223) 

(224) 

(225) 


‘dc . (l~k)df 

— (,o,r)~ -C2—0 


dr rf^i A,V 3 Jo(A„) 


64^ [l 


ia. 


4Ji( A n r) Anh-^f) , 

+ ^452 + ~ ~ ^ 

(226) 
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where 


64 ~ z^ 1 “ k ) c 2’ (227) 

Each of (222) and (226) can readily be integrated with respect to r since the integral of 

MKr) is -J 0 (A„r)/A n . This readily gives expression for the interfacial shape /(r) and 
concentration in the solid at the interface c s (z 0 ,r). 

We now turn to finding simplified expressions for the velocities near the R~ 1/6 boundary 
layer at z = z 0 for the simplified case of this subsection. Using (214), (215) and simplified 

expression for we obtain that each mode, i.e. coefficient of Ji(X n r) for the horizontal 
fluid velocity u scales as 


* R l/3 

p~^n(zp— Zj) r rjiO / \ i r-f ^'2 1 

tf’MXJTV 't«iX (*/) - ft.l [1 - r ). 


(228) 


The coefficient of J 0 (X n r) in such a series representation for the vertical fluid velocity w + 1 
scales as 


k /? 1//6 

Katt T p -*n(zo-z r ] 

x 1 /3 t o> 


- ft.] [1 - f ]• 

A n 


(229) 


7 Discussion of Results and Conclusion 


This is a problem with many non-dimensional parameters, as seen in Table 1. It is interesting 
to note that as R T - 00 , the leading order asymptotic behavior for the interfacial slope 
and radial segregation contain far fewer parameters, as seen in (212) and (213). However in 

applying this formulae, it should be pointed out that the next order correction in the result 
is 0(R C / Rt} or 0(R t ^ which ever is larger. 

To the leading order, we find that the interfacial slope and radial segregation scale as 


R 


, when other parameters are held fixed. Thus, in qualitative agreement with numerical 
results of Chang k Brown 16 , we find that moderately large gravity is worse than large gravity 
when it comes to controlling interfacial slope and radial segregation. 

To avoid a very detailed discussion of all the physical parameter dependence through m,, 
m 2 and &, we now restrict ourselves to the special case in subsection 6.1. It appears that 
most experiments satisfy the additional restriction placed in subsection 6.1 so there is not 
much loss of generality in doing so. 

The first notable observation is the crucial dependence of interfacial properties on the 
quantity z 0 - 2/ . This is so in the general case, where (212) and (213) hold, as well as in 
the special case where (222), (226), (228) and (229) hold. Physically, 2 0 - 2/ is the ratio of 
the distance of the interface from the end of the insulation zone in the solid to the radius 
of the cylinder. Since the smallest A n is A, = 3.83, it is clear that an arrangement that 


29 


would make * - a, (note here is determined from (37) when « « «/*> * 
large number, say 2 , will ensure that » at best «~, winch is rather s„ alh 

Since A, = 7.016 and other A„ even larger, it is clear that under most such , e^xpenm 
conditions, the dominant contribution to the series (214), (215), (222) an ( ) wi come 

from the » = 1 term. The profile of each of the interface /( r) and concentrate of solute 
c (r) will each look like .«A,r). In this context, it is mterestrng 10 " ot< - 1 a we ° ul ' 
that the numerically computed Chang & Brown- interface (or the Urge- « 

number quoted in their paper is in good agreement with a Jo( i>) P r0 seven 
numerical calculations deal with the fully nonlinear flow. The coefficients however could no 
be matched as their boundary conditions at r = 1 is different from ours. Company with 
the Adornato-Brown 17 profiles is not as favorable, especially near r - ■, presun 
the absence of an ampoule in our analysis. Returmng to relatrons (222), 26), ■ «• » 

that if most of the contribution in the series comes from the first term, ere w. e 

approximate proportionality between the radial segration & (r, so) and the mterfacal slope 
f'(r) with the coefficient of proportionality equal to 


(1-fc) 

D 


c 2 


1 - 


D {jd + \Jws + ^ o)j 


(230) 


Earlier, Coriell k Sekerka 22 hypothesized such a relation based on a purely diffusive “ lc “ la 
tion Brattkus k Davis 2 ', without any insulation zone, find no necessary relation. We fin 
here that with a proper insulation zone thickness, there is an approximate proport, ona] >ty 
between the two though the constant of proportionality in (230) ,s different from Gone I * 
Sererka 22 bv the appearance of the second term within the square parentheses. Since A, 
fairly large, if the Peclet number * is not large, then it can be expected that in an approxi- 
mate way, the second term within the square parentheses ,n (230) can be ^ >1 *“ d ’ 

which can again be small for many materials. This may explain why Cor, ell k Sererka go 

reasonable agreement with some experiment data. 

Now, let us discuss the physical meaning of the scale parameter * that appears 

(222) and (226). First note from (223) that <5, becomes large when 


T°‘ - (7 - 1 )mc 2 /D 
k 


(231) 


n ,1 i <5 also nresent in the denominator of £3 in (225). 

approaches zero. However, the term (231) is also presem m 

Tims when the expression (231) approaches zero, both (222) and (226) approach a finite 
limit. However, the expression (231) has to be kept positive to avoid constitutional super- 
cooling. Indeed, in the absence of capillarity, the condition for the onset of Mu hiis-Sererka 
instability of the basic r = 0 state is that the expression in (231) is zero. From the boundary 
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condition (13b), the quantity e[q lt T^( Zl ) - q 2 J is the negative of the radial temperature 
gradient in the solid at r = 1 , z = z/. Thus, for a dilute alloy, the ratio 


T 0 ' - - 1 )mc 2 /D 


(232) 


is approximately the ratio of the horizontal temperature gradients in the solid at r = 1. 
2 = z i to the vertical temperature gradient in the melt. Thus, h x = r t (« 5 / «)fl~ 1/6 , where 
k s /k is the ratio of the thermal difFusi vities between the solid and the melt. 

Now consider the parameter S 2 . It is clearly the 2 derivative of the log of heat loss rate 
at r = 1 , z = z/. If dimensional coordinate z were used, it would be the product of cylinder 
radius a and the logarithmic derivative with respect to z of the heat loss rate evaluated 
at the edge of the insulation zone in the solid side. This is clearly a property that can be 
controlled by appropriate design of the heat transfer properties of the Bridgman apparatus. 
When most of the contribution of the series in each of (222), (226), (228) and (229), comes 

from the n = 1 term, as is the case with an appropriately large value z 0 -z 7 , then the special 
choice 


S 2 = X x = 3.83 (233) 

will have the effect of minimizing the interfacial slope, radial segregation as well as fluid 
velocities in the boundary layer near the interface. When 8 2 < 3.83, the interface slope in 
(222) is positive, meaning that the interface will then be bulged towards the solid as reported 
in previous numerical computations (Chang & Brown 16 for instance). 

The role of the parameter S 3 as defined in (225), is more complex. For sufficiently 
dilute alloys, clearly since c 2 is small, 8 3 will be small. In that case, the term within the 
curly parentheses in (222) reduces to unity and therefore in that case there is no explicit 
dependence of the interfacial shape on the segregation coefficient k or the Peclet number 
D~ l except through 8 X and z 0 — z 7 . 

Even in a general case, the explicit Peclet number 1/D dependence in (222) and (226) is 
weak, except when Peclet number is comparable or larger than A] = 3.83 

We now discuss fluid motion. As discussed earlier, vigorous fluid motion is confined to 
the boundary layers near z = z 0 , z = z 77 and z = z 2 as well as near the side walls r = 1 for 
^ > z//. Elsewhere, the motion is 0(1) in the bulk for z > z 77 and o(l) for z < z 77 
as gravity is increased. This qualitative feature appears to be in agreement with previous 
numerical work. 

We now discuss details of the motion within the R T boundary layer of the interface. 
Here the scale of the horizontal and vertical fluid velocities for the special case of section 
6.1 are given by (228) and (229). Aside from the factor of (1 - S 2 / A n ), whose significance 
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has been already discussed, the remaining part of the scale factor, when written in more 
common notation, as explained in Table 2 is given by 

— R^ 3 Re- 1 Pr~ x r„ e - A " ( * 0-Z/) ( 234 ) 

for the horizontal velocity u. For the vertical fluid velocity tn + 1, the scale factor is: 

— R\! 6 Re- 1 Pr- 1 r s e -^ z °- zi) ( 235 ) 

K 

where * is the ratio of solid to melt diffusivities, fie is the Reynolds number, Pr is the 
Prandtl number and the parameter r s is given by 

efoi, 7?(*f) ~ 92m] (236) 

r s rpQ* 

Physically, r, is the ratio of the horizontal temperature gradient in the solid at r = 1, * = 
to the vertical temperature gradient in the melt corresponding to the Tiller et al solution 
Note r, differs from r, due to an additional term in the denominator. The requirement (220) 
for dropping nonlinear terms in the Navier-Stokes equation can be written as 

— jPr -1 r s e _Al ^ 0-27 ^ << 1 (237) 

k T 

However, we are of the belief that (237) is too stringent since within the boundary layer, 
the extra * derivative in the Stokes operator part of the Navier-Stokes equation introduces 
a factor that scales as fif which only needs to be larger than the convective terms, whose 
relative size is given by the left hand side of (237). This will have to be confirmed by a 
nonlinear analysis in the future. Nonetheless, it seems that through a safe choice of (z„ 
the requirement (236) could perhaps be met in experiment. Of course, z„ - raI,[ ''' 1 ® 

chosen arbitrarily large since large temperature gradients have to be maintained to avoid 


supercooling. 

Now, consider fluid motion within the boundary layer around 2 
horizontal velocity in (218) can also be written as 


zjj t The scale of the 


R]/ 4 Re- 1 Pr- 1 r, 


(238) 


where . cfar°(*/f)-ft) (239) 

7 rpQf 

is the ratio of the horizontal temperature gradient at r = 1, 2 = z// to the Tiller et al 
vertical temperature gradient in the melt. The vertical fluid velocities in (219) is smaller 
than the horizontal in this bounday layer by a factor of J# 4 , except near r = 1, where the 
fluid velocity components are of comparable magnitude. The requirement (221) for dropping 
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nonlinear terms may be harder to satisfy in an experiment, even though there is a smaller 
power of R t in the scaling compared to (237). This is so because there is no quenching factor 
e -Mz 0 -„) However, because of the boundary layer structure of the solution, the interaction 
of the interfacial properties at z = z 0 with fluid motion near z = zu or z = z 2 is weak. This 
suggests that one can have a highly nonlinear flow in the bulk of the melt, yet if the fluid 
Reynold number near z = z 0 is not large, i.e. (237) holds, then the asymptotic results (222) 
and (226) will also hold. In that case, it would be interesting to see if the condition (233) 
can be implemented through appropriate design of the Bridgman apparatus. 

Again, because of the boundary layer structure (i.e. exponential decay of the modes 
away from the boundary) of this problem, as it would be true for any thermally or solu- 
tally stabilized configuration, introduction of nonlinearity, when important, is expected to 
be mathematically manageable since by scaling the z - z 0 by R and using the scale infor- 
mation on velocities, a nonlinear boundary layer equation with only a few parameters can be 
formulated. Through appropriate matching with the linear solution in the nearly stagnant 
core of the insulated zone, a global solution can presumably be constructed. This will be 
subject of further investigation. 
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Appendix A: Derivation of effective interfacial boundary conditions 

The purpose of this Appendix is to to derive effective boundary conditions on the melt 
variables at a = * by using (105), (128), (118) and (119) together wrth the mterfac.a 

conditions (120)-(124). . . , / 19 o\ tW. for 

First, recalling that J? is a linear function as in (35), ,t « clear from (128) that 


° <Z<ZI ' am(z) = anr(z) + + Bre-'"* 1 -"', (j41) 

u„,M -Ih ..!? + *.]• ^ 

In the insulated zone 2 / < a < z„, on the other hand, 

a,„(z) = <43) 

The continuity of a„, and its first derivative at 2 = *, implies that each of B, and ft : can 
be expressed in terms of B, and ft. Further, from (118), B 2 can be solved rn terms of ft. 
Going through the algebra, we find the following express.on for «„ m the rnsu a e 


where 


< Z < Zq\ 


, a _ / >„(*-*/) _ - 2 \ n z, e -\n(z-zi)\ Bi+Gie^ Z,) +G 2 e An( * z,) a *np(0)e ; e 

— \ e ' 


— ^n 2 J a> ^n( 2 2 / ) 


where 


1 f , C , a 'snp( Z l) 

G 1 — ~ Q'snpyZl) ' \ 

Z A n 

If , v <np( Z l) 

Gi — n flsnp( 2 /) \ 


Using (A4) for a sn 


for a,Jz 0 ) and a'Jzo) and eliminating B, between them, we obtain 


(l' sn {z 0 ) — X\CL S n{ z o) 


where 


A', - >n (1 __ 

,j-G 2 + u„ p (0)e- A ""-G,e'“”''} 

-A„( 2 o- 2 /)i rr— • 


(1 4- e" 


V", = 2A n e 


1 — e" 


Similarly, the linear homogeneous e„uatio„ (105) with condition (119) can her solved 1m terms 
of one arbitrary constant. By eliminating this constant between 6„(so) and M-o). 

(,410) 
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where 



where 

V P3 — P4e Z °( pi 

y\ 2 — * 

1 _ e*o{P4-p3) ’ 

(All) 


P3 2 D s 1 ]/ 4 D 2 S 1 A "’ 

(A12) 


Pi 2D, + A “' 

(A13) 

From (120)-( 1 22) (with d 0 = 0), we can express 



_ «"(*>) + TjM*>) a n (z 0 ) + TJKizo) 

T s' + mc° s '(z 0 ) T 0 ' + mc°'(z 0 ) ' 

M 14) 

Equation 

(120) then becomes 



a n{ z o) = ^11<2sn(-2o) + n l2^sn{ z o)i 

(A15) 

where 

T o- _ JO' 

nil = 1 2 

zy+? <w 

(A16) 

ma T 0 ' - T°'l 

«12 = f L 

W [T? + fc°'(z 0 ) J 

From (122) and (A 14), we get 

(A17) 


^n(^o) — ^21^5n(^o) + ^22^n(^o) 5 

(A18) 

where 

„ _ fi [<*'(*,) - lcf( 2o )] 

“ |I? + f <?(*>)] ’ 

(A19) 

Thus, 

n22 = I + m KM - KM 

k i |r» + =< ( 2o )] 

(A20) 


a,„( 2 „) - ’ h2a '‘(~°) - n ut>n(zo) 

n 22 n n — n 12 n 21 

(A21) 

A ( „ \ n 21^n(^o) + n\\b n {Zo) 

°sn{Zo) ■. (A22) 

n 2 2 n n ~ n 12 n 2 , y ’ 

Further, from the simplification of (123) into <(*,) = £<(*) (due to linearity of 7? and 

T°) and (A7) and (A21), it follows that 


n f o) + Mn + /? 2 &n = 

A n 

(A23) 
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where 




XiK,U23 

A n K(n 2 2 n 11 — n 12 n 2 l) 


03 — 


AA K 5 T 212 

^ 2 A n /c(n 22 nn — «i 2 « 2 i) 

2A„/c s [ — G 2 4- a»np(0) e ~ AnZJ ~ (ji e ~ 2AnZf ] 
ac [1 - e _2A " Zo ] 


(A24) 

(A25) 

(,426) 


Using expressions for n,„ n,„ »„, n 2J , X„ G, and ft in the above and carrying out the 
algebra, we get m c°'U))(l + «-»-») (427) 

Ac(mc»'(2o) + T°')( 1 - e _2An *°) ’ 

Ac a mct(T a °' — T°')(l + e~ 2AnZ °) (A28) 

^ 2 ' _ n(3(mc°'(zo) + T °‘) (1 - e" 2 ^) 

+ e- 2A " z 0 ~ 2 (qWT^zi) ~ 92,)e- A ^ - + e^ 1 )} 

K [1 _ e -2A„zo] (,429) 

ing (124), (A14), (A22), we obtain another effective boundary condition for the melt 


P 3 — 


K s {(quT?(zi)-q2,)^ 


Using 
variables of the form 


where 


b' n (z 0 ) + m 1 6 n (zo) + m 2 a n {z 0 ) - 0, 

mi = 04 + ^205* 
m 2 = 06 + -^207 


In the above, 


* 

a = + 


$3 - 


, m (%rf (zo) + 1 Tr 1 n 0 '(*o) - n 0,, (ro)) 

(T 0 ' + mc°'(zo)j 
D,nn 

~ D(nnn 22 — «i2«2i) 

0 (%-Ca"(gp) + ( Z o) — C ° ( Z °)) 


(A30) 

(-431) 

(A32) 

(A33) 

(434) 

(A35) 
(A 36) 


a (T 0 ' + mc°'( z o)) 

£>»na i 

^ 7 D(nn n 22 — «i2«2i) 

r ,nrl « in the above the expressions for 05 and 07 simplify 

Using expressions for nn, rin, n 2i anc * n 22 ’ 

t0 ' D a (mc°'(z 0 ) + kT°) (A37) 

= D(m^(2 0 ) + r° , r’ 
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g = ££ q) - <?(*)) 

Da(mc°'(zo) + T 0 ') ’ (^438) 

Equations (A23) and (A30) form the effective boundary conditions on the melt variables 
that incorporates all the the coupling between the variables on the solid and melt side. 
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Appendix B: Derivation of the asymptotic form of approximate general solut.on 

The purpose of this Appendix is to find approximate expression for c„„ a,,,, K„, F i 

and Ft appearing in (146)-(152) Rr » 1 with l«d « Rt ' 

ft is possible to find uniformly valid asymptotic representation of the S ° " fc 

this limit. Nonetheless, tills expression is too complicated and not very sui a e 

• if \ n R << R b/ 6 So we first derive easier expressions for the asymptotic. 

manipulation when R c « R r • (\*R T f. Clearly, even when 

behavior of the solution for R T » 1 and K\ R c\ \ n ) 

| | « lif\ this condition is not uniformly valid for all A„. Thus, for arge enoug i », 

L have to use the more general expression-, nonetheless, from the arguments ... section 6 
we know that the detailed behavior of the solution for large A„ « not necessary to conclu 

that the contribution to the series (96)-(101) from large A„ » mdeed negligible and so 

, , ii- for I ft 1 < 6 can be deduced from the results in the following 

behavior of the solution tor \R C \ < cau uc 

subsection 

Bl. Subcase A*|/2 C | << 7 5 > R t >> ^ 

From (143), in the heat zone 2// < 2 < 20 , we have 


[£l£ 2 £ 3 + >< 2 n R cC\h(z) - XlRrCl ] C » - ^2 [ A »( q\T + flfe)] • 

To the leading order, we expect 


C n ~ C n> 


[£,£,£, - XIRtCi ] c„ = £2 [-A n (-4,r° + 92 )] ■ 


(Bl) 


(B2) 


(S3) 


where c° satisfies 


Thus, we can write 

[£,£ 3 - A iRr ] c« = -A n (-„T 0 + 92] + El -2 {i-lfiliab,] - V" ( ' ^ (B4) 

where A\, A 2 are arbitrary constants and 

7 = A )/ 3 4 /6 (B5) 

J r l 2 \i (^ 6 ) 

U [y] = y - A « 

Li[y] = D[y 2 - A nl + y 

l 3 { y] = ( y 2 - A ') 2 + l(y 2 - (B8) 

1 ,/jlTa* ^ B9) 

^'2Dl4D !+n 


40 



On solving (B4), we obtain, 


P2 2£> + V4Z) 2+A ” 


where 


c n = c°np + 52 7 2 A j e p > (z ~ z °'> + £ Aje^-to), 
f =1 j = 3 


„0 _ A n [ gl? 10 + <72] X h n q\T 0 ' 


np [A^r + A®] v[\IR t + A«P’ 
and 7 j are the six independent roots of 

MTiJMTi] - A^r = 0 

that are labelled such that for R T » A*, 

7j ~ 7 w i 

where w,- are the six roots of unity defined by 


u>3 = —e ,ir / 3 = —uj 8 


U) 4 = — 1 = — (jj 7 

U> 5 = — e _,1r/3 = -U7 6 


In order to find the next order correction in the asymptotic 


expansion 


c n ~ 


C « + Cl + ... 


we notice that c\ satisfies 

A [AA - A*/? T ] 4 = -A^ c £,[/icy-A^ c7 - 2 

where 




w=R -p 


7j 7> ■ 


Then one can solve (B19) in a standard manner to find that 


4 = 4,-7- 2 a 


'4jTi[pj]/(r:)ePdr-ro) 


(BIO) 

(B II) 
(S12) 
(B13) 

(B 14) 

(Bli) 

(BIO) 

(Bn) 

(BIS) 

i= 3 

(B19) 

(B 20) 
(B21) 
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„ * (B22) 

^3 Miil - A » fiT > 

where c t \ p is a particular solution of the form 

c' np = ( A 2 {z - 20 ) + AOfe(z) 


that satisfies 

C 2 [-Cl £3 - X^Rt] C np = -^n^c£l[* C np]- 
When >> A*, the above asymptotes to 


n r>*j 

C np 


9< 2 ) 


B ( , 32 c 

« f t,[-F i{ tll ' r| ‘ W D C ' ,+ LA-D-'\\ 


Thus, if we write 


c„ = c„ v + j:r‘‘C,F' + 'tc,F' 

j = 1 J = 3 


(5 23 ) 


( 524 ) 


as in (146), then 

Cnp ~ C np T C np + 

. -i p \4 r an d its derivatives are (^(/Ix 1 ), which follows 

Note that m the special case K T « A „, c n P ana «* ucl1 / 

from (B12) and (B23). This scaling property is used in section 6. Further, in ( 

f J, , \ f™- 7 = i9. 1525) 


pc ~ e nr,(*-*o)|l 
pc „ ^,(*-22) 1 1 _ 


\ 2 n RMlj}h(z) 


L 2 [l j] {L\{l]}L 3 {^, 

XlRcLSMz) 


M z ) | \ 

Til -XIRt) "J 

-1 


for j = 3,4,5 , (526) 

-f..l for j = 6,7,8 , (527) 


Z. 2 [Tj] [M7iM7j] - A n^r] 

where we choose new arbitrary constants 

C, = 7 s . 4 ,, C 2 = 7 s A 2 e« l "-“>, Cj = 7 5 A, for j = 3 , 4,5 C; = 7 5 Aj< 4 ' ,< »— > for J = 6 , 7 ^ 

Using (145a), we can solve for a„ once c„ is known. Such a solution also needs to be consistent 
with (144a). We find that 


where 


a„ ~ a np + ^7 l CjFj , 
j=i 


o i „1 i 

anp ~ a „p + ®np + — ’ 


( 529 ) 

( 530 ) 
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where 


„o _ ^[-< 7l 7’° + <?2 ] \l qx T 0 ' R t 

np KRt + a« v\KRt + KY ’ 

and a] ip is a particular solution of the form 

< = \B*{z ~ Z 0 ) + A] h(z) 

that satisfies 

a lp = ~^ 2 n Rrc\ p ■ 

When R T » K ’ the expression simplifies so that 


,1 


XlRch(z) 


Qnp ~ l2hD-']Li[-D-'] 


i'ri-D-'K + 


Lrl-D-'} 


np L 2 [-D~'] np 


{B 31) 

(532) 

(533) 

(534) 


Further, 

F a e Pj( z ~ z { 2j-2)) 


F^ ^ * — p^j( z z o) 

3 ^2 _ \ 2 e 

ij 


, gcg/MfjMgj] (-I 2 fp,-I5 3 [pj] + 5A^5 t ) \ 

DR T p j L 2 [p j ] (h \pi\L 3 \pj] - XI Rt) + ") ° r 1 ’ 2 ’ 

(535) 

_ MjMMM + DA »^) 


55T7jL 2 [7j] (Ii[7>]Z/3[7j] - A*5 T ) 


+ .. for j = 3,4,5 , 


TTta 

r • ~ 

3 ^2 


^ e'Yjf-e— ^ 2 ) 


2 _ izMMMM + ^ a ^t) \ 

55t7j ^ 2 [t>] (5 1 [7/] r 3 [7i] - 5 r ) + "/ 


i]-xi 

Using b n = —C^Cn — a n , we can write 

bn = 6 n , + £ 7 4 C,P‘ + A iKT'i.CjF* , 

where 


i=3 


(536) 

+ .. ] for; =6,7,8 . 

(537) 

(538) 


^np ^[Cnp] (Inp ■ (539) 

Simplified expressions are possible for R T » A, 4 ,, which we do not care to write. The only 
property that will be of importance in the analysis in section 6 is that b np is this range of A n 
is 0(R c /R t , 5f*), as is a np , which follows from (B30), (B31) and (B34). 

Further, in the general case in (B38), 

5 fc ~ ( Pi L z\Pj] , Rch(z)L x \p 3 ] Pj \ . . . 

I Dy‘ ~R^ TMh 1 ,= ' 2 ' (S40) 




(B41) 
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(542) 


n 


_ ( =^j( 2_2 2 


-«)l2hMh(z) + .. for j = 6,7,8 

* r f * 1 ' ' 


Now, we comment on the validity of the asymptotic behavior given so far. The behavior 
of the particular solutions K r and is clearly valid anytime \R C \ « Rt Tins is 
uniformly true for all A,,. Further, expressions (B25), (B35) and (B40) is also consistent m 
this regime for all A„. However, the expressions (B26) and (B27) do not remain valid when 
lAjfcl is the same order or larger than y 5 , i.e. (A IRtY'*- This is because from defimtion, 

L\bj]Lzbj] = x I r t = 7 6 - 

Thus, for large Rt, 

L\bi\ L *bi\ = °b*)- 

In order that the second term in (B26) and (B27) be smaller than the first term, it is necessary 
that | A 2 R c \ « 7 5 - This is the origin of the restriction for the asymptotic behavior shown 
thus far. Note that the condition A>|ft t | « T cannot be uniformly valid for all A„ even 

when \R C \ << R T ■ 

B2 General form of Solution for |5 C | << Rt 

In order to find uniformly valid expression for all A„ and at the same time find expressions 

valid for \RA « Rt, we find six independent solutions to the associated 

equation in (Bl) in the WKB form. These solutions will replace the expressions ( I an 

(B27). Once this is found appropriate expressions can easily be found to replace 3 ), 

(B41) and (B42), which are also invalid in the general case. 

We consider WKB solution of the form. 


P'9 = e W} for j = 3,4..,8 


(543) 


We will think of the relation (B43) as exact and relate the WKB approximate behavior 
through the relation daa\ 

Wj ~ w 0j + w h 

We find that for fir » 1, with \Rt - R<\ » 1 , a uniformly valid expression for is 

given by the following expressions 

R 


W 03 = ~ J*[{ X l R T - X ljf h ( Z )) Lie ' 23 + X nb dZ 

Wo 4 = - J\( x l R T - x l^H z )) 3 + x l\ kdz 

W 0s = - [ Z {i X l R T - X ljfH Z ))' Je ~ t2 * + X nb dz 

J Zo *** 


(545) 

(546) 

(547) 
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w o* = jjXlRr - A l*kh(z))$e~ i2 S + A 2 Jdz 

^°7 = ~ ^n~jjh(z))z + A 2 n ]*dz 

Rr 


(5 48) 
(549) 


W ° '» “ J Z2 - \ 2 n ^h{z))h' 2 l + \ 2 n ]Uz ( 550 ) 

In the above, we proceed with the understanding that the principal argument is being used 
in taking the squareroots and cube roots. W 1} for 3 < j < 5 is determined by 

W h = f o W(.(z)dz ( 551 ) 

and for 6 < j < 8 , is determined by 

W h = (552) 

where each W" are determined in terms of the corresponding through 

= w « { 27 « - A n) 3 + 24 a;« - a;> 2 } + (§ + i) - a?j 3 - “i 




(553) 


- K) 

Further, note that when R T » |5 C | with A* « R T , these expressions reduce to 

W o } ~ (556) 

in agreement with (B14). Thus, the expression (B43) is a generalization of the earlier 
formulae for F/ for j = 3, 4, ..8 

Now consider, finding F/ for j = 3,4..8. It is determined as a particular solution to 


C X F° = -\ 2 n R T F< 

Through a standard dominant balance procedure, we obtain for j = 3, .. 8 , 

1 


ZT'a 6 Wi 

r ~ — 7 e J 


W’ 1 


,w ?-\ 2 (Wf-xiy 

Since Fj = -F/ - £ 3 F/, it follows that for j = 3,4, ..8, 

-1 


+ 


-1 


(557) 


(558) 


F 6 ~ 

J Wf - A l 


((Wf - A 2 ) 3 - 7 6 - 2(Wf - A 2 ) 2 !*;" + 4IF 7 'X /, (^ ~ *1) + ^(Wf - A 2 ) 


(559) 
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Notice that it is necessary to retain all the terms in (B59) even when the first two terms 
within the parentheses are clearly 0 < 7 6 ) while the remaining are 0(7 ) because there 
cancellation between there is cancellation these leading order terms because 

(< - ^) 3 - 7 s = 

and this need not be larger than -f. Relations (B43) with asymptotic behavior (B44) and 
relations (B45)-(B53) replace the more restricted express, on (B25) and (B26) wen 

0( 7 5 ) or larger, provided Rt » W- The remammg expressions (B25) (B12) re, 

valid. As far as «„ and 6„, expressions (B57) and (B58) replace (B36)-(B37 and (B4 )•( ) 

respectively, other expressions (B35), (B30), (B31), (B32), (B39) and (B40) st.ll reman, 

as long as Rt » |/u 
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Appendix C: Asymptotic evaluation of solution to MZ = S 

The purpose of this Appendix is to carry out the asymptotic evaluation of solution to 
(184) when R T » 1 and A n = 0 ( 1 ). In this case, 7 defined by (147) is » 1 and the 
asymptotic relation (B159) for the roots 7 j hold. 

First, we consider the elements Pj through P 8 that appear in the expression for S in 
(186). Since for 1 < l < 4, the elements A/ u e -lv *^//) sca j e ^ e -y 3 {z n -z 0 ) for j _ ^ 7 ^ 

and is therefore transcendentally small in R T . We ignore such terms in P h Similarly, 

for 5 < / < 8 , the elements M Uk e~ w ^ z ") scale as e^~ z ") for j = 3 , 4,5 and is 

transcendentally small in R T . This leaves us with the contribution from terms multiplying 

Hl ' J and ln ( 187 ) and ( 188 )- Thus > it is necessary to calculate the first two rows of the 
matrix //, which is the inverse of G. 

First, consider the simplification of the matrix G, whose elements are defined by (175), 
(176), (179) and (180). On examination of (153)-(170), it is clear that 


G ~ G° + O (\lR cl ~>) , 


where for k = 1 , 2 , 


„k— 1 


(Cl) 


(C 2 ) 


Pj ^ j - 1, 2, G£j = — P for j = 3, ..8 , 

7j ~ -^n 

while for k = 3, ..8, 

«(* -2 ) / \ (k— 3) 

k ’ j = Z) 7 (^-2) f or J = 1)2, G^j = ^ j for j = 3, ..8. 

Now consider the problem of determining the the first two rows of H° = G°~ l . This 

be conveniently done by finding the first two components X[ and X l 2 of the vector X 1 that 
satisfies 


(C3) 


can 


G°*' = £' = ( 0 , 0 , 0 ,,, 1 ,.. 0 )\ 


(C4) 


the only nonzero element of E‘ is a one at the l-th entry. It is convenient to define symbols 
r k so that 


Thus, 


isM = (y 2 - Aj) ! + V-A;) = £ rt / 

k=0 


r ° - A » ’ r » = —f » . r 2 = - 2 A, 2 , , r 3 = - , r 4 = 1 . 

V V 

Also, we define r h = 0 for negative k. Notice that (158) implies that 

X 2 = 7~%[7j] 

h n 


(C 5) 


(C6) 


47 


So, when the first two components of the vector identity (C4) is written long hand, we get 


4 8 


x[ + x‘ 2 = E i~ 5L M xt i + ^ 5 £ r * S ^ 

j=3 *==° t=3 

p,*j + p 2 X< = E 7 " 5 7i = 7' 5 E r* E 7,<* + 0*{ + > ( C8 ) 

j_3 fc=0 J-3 

where $ fcj - is the usual Kronecker delta symbol. From the (k + 3)-rd element of the vector 
identity, with k ranging from 0 to 5, it follows that 


E7‘*'=-^ep‘ +i *;+**•« 

j=3 l U j - 1 


Using this in (C7) and (C9), we get 


X[ + X[ = -TPS E Pj L M X j + 7 &r '-3 + *u> 
^7 ;=1 

Pi *{ + p 2 x' = + ^' 5r '- 4 + *'- 2 ‘ 


from which it follows that 


hi, = x! = [i + ^.w] 


, (P2n~3 - n- 4 ) 

(P2 “Pi) 


H\ l, = A= 1+^iilPi 


Pi Stj — hg _j_ _s (pi r <-3 ~ n- 4 ) 

Pi — p 2 (pi — P 2 ) 


(CIO) 


(Cll) 


(Cl 2) 


(Cl 3) 


Now, let’s consider the other elements //°j for k > 3. Clearly, from (C9), it follows that 
that for k ranging from 0 to 5, 


E 7 _fc 7j Hj,t = — + «W-3 

j=3 J =1 


(C14) 


Using the results (C12) and (C13), it is clear that the right hand side of (C14) is 0(7" 2 ) for 
l _ i 2 and 0(1) for / = 3,4, ...8. Since the coefficient of each term on the left hand side of 
(C14) is 0(1) in the asymptotic limit 7 -* », it follows that 

H] t = 0{ 7" 2 ) for j = 1,2 and H?, = 0(1) for j = 3, ...8 (015) 

Now, since from (174), (B12), (B30), (B31), (B34), a np and its first derivatives are 
0(RJRt , ftr’ )’ while C "P ancl its derivatives are 0(7^*), it follows from (174) and (178) that 
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each of and R 2 are 0(R c /R T , 7 6 ), while R 3 and R 4 are 0 { 7 -1 ), and 0(7“ 2 ) respectively. 
All other components are 0( j~ 3 or smaller. This means that for j = 1,2, 


le -wA*n) £ H hk R k = x 0 ( 7 - 5 


MR c ,r 5 


). 


(C 16 ) 


and so 

P l = e - ^-^) X 0( 7 - 5 A^ c , 7 - 5 ) for /= 1,2, 3,4 . ((717) 

P, = c pi (*»—•*//) x 0 ( 7 -5 a 2 7? o7 - 5 ) for / = 5, 6, 7, 8 . ((718) 

Further, we notice that in (186) in the terms beside P t , all except S 3 are 0( 7 _5 /7 c , 7 -1 ) = 
o(l), where as S 3 contains a term that scales as which will dominate every other 

term for A n = 0(1) and z 0 -z,= 0(1). Thus, to the leading order 


5 ~ S° = A;% e - A "^-*')(0,0,l,0,0,0,0) (019) 

Here we comment that even if we were not to ignore the contribution to Si for / > 5 
on the grounds that z 0 - zj is large and therefore small, it would not affect our 

leading order results for Zj, Z 3 through Z 5 because of the the special structure of the reduced 
matrix M in this limit. Therefore, the results quoted in (209) and (210) would equally be 
valid in this case since the contribution of Z 2 , Z 6 through Z 8 for interfacial properties a n (z 0 ) 
and b n (z 0 ) (given by (199) and (200)) are exponentially small. 

Now, consider simplification of the matrix M, whose elements are shown in (189)-(198). 
We notice that 

M ~ A/° + 0( 7 -1 , 7 - 5 A 2 fl c ) , (C20) 

where the only nonzero elements of M° are 


1 for k = 1, 2 j = 3,4,5 , (<721) 

M h = , ((722) 

= _Pl “ m ! + m 2 , A / 4 °2 = e P2(zo ~ z ^ (-p 2 - mj + m 2 ) , = -m 2 u)J 2 for j = 3, 4, 5 , 

((723) 

= ePl( "“ 20) , <2 = 1 , Ki = ~“7 2 ** j = 6,7,8 , ((724) 

M 6J = 1 and M 7 j = u>j for j = 6, 7, 8 , ((725) 

<1 = -e Pl( ^o) , < 2 = -l. ((726) 

The solution to 

A/°Z° = 5° ((727) 
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can 


be found in closed form. We find (Z%,Z%,Z%) — (0,0,0), 

„ / 1 -v/3 , 1 V3 N 

(ZlZ°,Z°)=[-- + ' 1 r, 


fhc 


-A n (z 0 -zj) 


Once this is found, it is clear that 


‘Im-ifae 


— A n (-0 z l ) 


(Z°,Z°) - - + _ e ( Pl - P 2 )(z2-zo)(p 2 + m, - m 2 )) 


(l,-e pl(22_2o) ) • 


(C 28) 


(C29) 
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Figure and Table Captions 


Figure 1: Schematic diagram of the Bridgman apparatus. 

Table 1: Nomenclature. Note all listed variables are non-dimensionalized. 

Table 2: Relation of non-standard symbols with more common notation in Literature. 
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Table 1: Nomenclature 


c 

(u,w) 

z 

zj 

zu 

zo 

Z2 

r 

a 

P 

v 

K 

D 

D* 

^1. 

f<72 

^2. 

C2 

k 

V 

9 

St 

do 

m 

Rt 

Rc 

h(z) 

An 


zJjl where T: temp.; Tq: planar interfacial temp for pure material 


Nondimensional temperature, ^r_f o 

and 7 2 : top temperature, 
relative concentration of one binary component 
Fluid velocity relative to velocity of container, u and w : radial and axial components 
z-.h. f where z: axial coordinate and z\\ z at the lower cylinder end. 

Value of 2 at the lower end of the insulation zone in Fig. 1 
Value of z at the upper end of the insulation zone in Fig. 1 
Value of z at the unperturbed interface 
Value of z at the cylinder top in Fig. 1 
- , where f : radial coordinate. 

Nondimensional thermal expansion: a = a(f 2 - To), where a: volumetric thermal expansion coefficient 
Coefficient of volumetric expansion due to increase in c 
Actual Viscosity/(a[/„) , where a: cylinder radius, U a container velocity. 

Thermal diffusivity in the melt /( aU a ) 
solutal diffusivity in the melt /( aU a ) 

Thermal diffusivity in the solid /(at/ a ) 
solutal diffusivity in the solid /( aU„ ) 

Biot number related to horizontal heat transfer on the melt side 
, where t « 1 is chosen so that q\ = 0(1) 

Biot number related to horizontal heat transfer on the solid side 
, where e << 1 is chosen so that q\, = 0(1) 

Nusselt number related to horizontal heat transfer on the melt side 
Nusselt number related to horizontal heat transfer on the solid side 
Concentration at the cylinder top 
Seggregation coefficient 

Fluid pressure w j lere p Q - density at temperature To and concentration c 2 
p oU a 1 • / r r 2 

acceleration due to gravity x a/U a 
Stefan number: 

surface tension/(Latent heat x a) 

Slope of liquidus line x (T 2 — To) 

Thermal Rayleigh Number 2 ^- 
z derivative of T°, the unperturbed T corresponding to c = 0 

Solutal Rayleigh Number jfc ) 

where c?; interfacial concentration of e = 0 state 
* — 

Convenient symbol for e n 

nth positive zero of Bessel function J\ 

Convenient symbol for [\„Rt] ^ 


where 


iO ; . 


Table 1: Nomenclature. Note all listed variables are non-dimensionalized. 
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Table 2: Relation of non-standard with standard parameters in 

Literature 


= Pe 1 , where Pe: Peclet number 

= (D s / D)Pe ~ l , where D s /D : Ratio of solid to melt solute diffusivity 
= where Re: Reynolds number based on translation rate U a - 

— i4“ J , where .4: Ratio between cylinder radius and length. 

= Re~ l Pr~ l , where Pr — Prandtl number 
= (k s / n)Re~ l Pr " 1 , where /c 3 /k: Ratio of solid and melt heat diffusivity. 

Thermal Rayleigh Number where~ denotes dimensional quantities, 

T 0 : Limiting linear Temp, profile of Tiller et al (1953) for /c/z 2 >> T 


Solutal Rayleigh number g ^ ^ 1 (4 - 7 D Jp 7 y > where c,: Tiller et al (1953) interface concentration. 

= flr Re~ 2 Pr~ x A -1 , where A = y ■ 

= ftTfte-’Fe-^c. -czJ-’A- 1 , A = 1 _ 


« 


Table 2: Relation of non-standard symbols with more common notation in Literature. 
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